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A convenient formalism for averaging the losses produced by gravitational radiation backreaction 
over one orbital period was developed in an earlier paper. In the present paper we generalize this 
formalism to include the case of a closed system composed from two bodies of comparable masses, 
one of them having the spin S. 

We employ the equations of motion given by Barker and O'Connell, where terms up to linear 
order in the spin (the spin-orbit interaction terms) are kept. To obtain the radiative losses up to 
terms linear in the spin, the equations of motion are taken to the same order. Then the magnitude 
L of the angular momentum L, the angle k subtended by S and L and the energy E are conserved. 
The analysis of the radial motion leads to a new parametrization of the orbit. 
ON . From the instantaneous gravitational radiation losses computed by Kidder the leading terms and 

the spin-orbit terms are taken. Following Apostolatos, Cutler, Sussman and Thome (ACST), the 
evolution of the vectors S and L in the momentary plane spanned by these vectors is separated from 
O |i the evolution of the plane in space. The radiation-induced change in the spin is smaller than the 

leading-order spin terms in the momentary angular momentum loss. This enables us to compute 
the averaged losses in the constants of motion E, L and Ls = Lcosk. In the latter, the radiative 
00 ■ spin loss terms average to zero. An alternative description using the orbital elements a, e and k is 

given. 

£N| . The finite mass effects contribute terms, comparable in magnitude, to the basic, test-particle spin 

^ ' terms in the averaged losses. 
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I. INTRODUCTION 



o. 

Coalescing binaries are important and copious sources of gravitational waves for the projected interferometric grav- 
CTN ■ itational wave detection experiments LIGO, VIRGO and LISA. The requisite signal templates challenge theoreticians 
I to predict the behaviour of such systems under gravitational radiation backreaction. Much progress has been made 
^J^ 1 in this direction mainly by perturbative approaches. In the initial epoch, the separation r of the two bodies forming 
the binary is large relative to the Schwarzschild radius and the relative motion can be considered slow. This enables 
one to employ the post-Newtonian expansion parameter e « v 2 « m/r Alternatively, with the system of units 

chosen appropriately, the inverse 1/c of the speed of light can be used as an expansion parameter JpQ. Computations 
of instantaneous radiation losses of energy, momentum and angular momentum up to (post) 5 ' -Newtonian order, 
including the spin-orbit and spin-spin interaction terms, were given by Kidder Q . 

The equations of motion of spinning bodies, disregarding radiation backreaction, were first considered by Barker and 
O'Connell Q, then by Thorne and Hartle ||. The relevant equations can be derived from a generalized Lagrangian 
depending on relative position, velocity and acceleration. The two-body problem is reduced to a one-body problem 
by eliminating the center of mass. The masses mi and mi still appear in the Lagrangian as reminders of the initial 
two-body character of the problem. Apostolatos, Cutler, Sussman and Thorne describe the evolution of the spin 
and orbital angular momentum vectors in the presence of radiation backreaction. They supplement the equations of 
motion with the leading-order averaged gravitational radiation loss terms, computed by Peters and Mathews [§,[§. 

In the test particle limit the leading spin terms from the averaged losses in energy, magnitude of the orbital angular 
momentum and spin projection of the orbital angular momentum, respectively has been computed by Ryan pi for 
generic orbits. Generalizing his earlier work on circular orbits |I~if| and Shibata's work on equatorial orbits E3 , Ryan 
obtains the first-order spin corrections to the Peters' equations of radiation-induced change in the orbit parameters. 
In a previous paper , to be referred to as I, we completed the description of the test particle by giving the averaged 
losses in terms of unambigous conserved quantities and by computing the radiation induced change of the remaining 
orbit elements. Our description relies on the use of Eulerian angle variables, a new parametrization of the orbit and 
averaging by the residue theorem. 

The computation of finite-mass effects can be carried out by choosing a suitable spin supplementary condition 
(SSC). Three convenient choices of SSC are discussed by Kidder A treatment of the nonradiative evolution using 
a noncovariant SSC (A2b) of j| can be found in Wex's work , and some radiative losses in Rieth and Schafer jl5) . 



1 



The finite mass effects in the spin terms of averaged losses, using the covariant SSC, were considered previously by 
Kidder, Will and Wiseman [j| and by Kidder [||. They computed the energy loss and the angular momentum loss 
for the case of circular orbits. It is the purpose of the present paper to obtain, for generic orbits, the averaged losses 
of energy, magnitude of the orbital angular momentum and spin projection of the orbital angular momentum. 

We employ the second order Lagrangian suitable for the description of the motion of comparable mass bodies in 
the presence of spin, following Ref. We want to describe finite-mass effects in binaries consisting of a compact 
body captured by another, massive, spinning body. The latter is exemplified by a spinning black hole or neutron star. 
In this way we generalize the picture of the evolution of a system consisting of a black hole and a test particle, which 
we considered in I. The equations of motion and the spin precession equations valid to the (post) 2 -Newtonian order 
are reviewed in the second section. The constants of motion are the energy E and the total angular momentum J. 
Since our goal is to get the radiation losses up to the leading terms in the spin, we approximately picture the spin 
precession by a rigid paralelogram with sides L and S, rotating about the total angular momentum vector J. In this 
approach the magnitude L and the spin projection L$ of the orbital angular momentum are constants of the motion. 
From their expressions an uncoupled radial equation can be found. 

In the third section, we determine the turning points. We then introduce two parameters, the 'eccentric anomaly' 
and the 'true anomaly'. They have the respective properties, and in the no-spin limit reduce to, the eccentric anomaly 
and true anomaly of Kepler orbits. The parametrization of the orbit by the eccentric anomaly enables one to compute 
the orbital period defined as twice the period of time elapsed between consecutive turning points. 

In the fourth section, we compute the instantaneous radiation losses. For this purpose, we make the assumption 
widely used in the literature that the radius R of the spinning body is comparable with its mass mi . Then S/r 2 w 
miRV/r 2 w e 2 V and L/r 2 w fiv/r w e 3 / 2 . Thus the spin is of e 1 / 2 order smaller than the orbital angular momentum. 
Under the assumption that the rotational velocity V is large, Apostolatos, Cutler, Sussman and Thorne H estimated 
the 'S-changing piece of the radiation-reaction torque' (the radiation loss of the spin) to be of order re R / 2 . This is 
smaller by two full post-Newtonian orders than the leading term in p6| dJ / dt which is of order re 9 / 2 , and by one half 
order than the spin-orbit terms given by Kidder 0j which are of order re 6 V. At the end of the section we present an 
estimate of the order of magnitude of the radiation loss of the spin for smaller values of V. We then find that, to the 
order we are interested in, we may interpret the instantaneous total angular momentum loss as consisting entirely of 
the loss in the orbital angular momentum. These considerations enable us to compute the loss of magnitude and spin 
projection of the orbital angular momentum. However the latter will still contain radiative spin loss terms, which are 
evaluated from the Burke-Thorne potential. Fortunately, introducing suitable new angular variables, all losses can be 
expressed in terms of the true anomaly parameter alone. 

In the fifth section we average these expressions by using the residue theorem. We find that the spin-loss terms 
give no averaged contribution to the losses. The finite-mass backreaction effects are of the same order as the other 
spin effects. The results obtained in the Lense-Thirring case emerge as smooth limits, although the test particle case 
is outside the framework of this paper. 

The losses are given both in terms of conserved quantities and in terms of 'ellipse parameters' a, e and the angle 
k between the spin vector S and the orbital angular momentum L. We give also the changes of the geometric orbit 
elements a, e and k due to radiation backreaction. This is carried out in the sixth section. 

As in ACST, we separate the motion of the spin and orbital angular momentum vectors in their momentary planes 
from the motion of the plane. Our approximation in describing the radiation effects on the former motion up to linear 
terms in the spin. 

Since we eliminate the post-Newtonian parameter v 2 from the losses at the very beginning of the averaging process, 
we keep the gravitational constant G and the speed of light c in our formulae. Our post-Newtonian parameter is 
e w v 2 /c 2 « Gm/c 2 r. The explicit c dependence allows an easier bookkeeping in the post-Newtonian expansions. 



II. THE ORBIT OF THE BINARY SYSTEM 



We consider the motion of a bound two-body system with masses mi and mi and spins S and 0, respectively. 
Considering only the leading-order spin-orbit coupling, and adopting the spin supplementary condition of pL the 
Lagrangian is: 

C=^ + ^+SC, (2.1) 
2 r 

where r = |r| is the relative distance, v is the relative velocity, and the perturbation term due to rotation effects is 

SC = 2(1 tf^ v(r x S) + ^v(a x S) . (2.2) 
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Here /i is the reduced mass and m the total mass of the system, 

m 1 m 2 m 2 
a = , to = toi + TO2 , r\ = — . (2-3) 

TOi + TO2 mi 

The Lagrangian depends on the relative acceleration a. The application of a generalized, second-order variational 
formalism yields the equations of motion: 

a = -^r + -£L J ^+^ r [ (r x v) • S] - (4 + 3r?)v x S + 3(2 + ?j)-r x si , (2.4) 



and the momenta: 



p.ae-t.^+e+^xs. (2.6) 

The energy i? and the total angular momentum J are constants of motion: 

B = p . T + q . a . £s f ^ + &. (rXT ), (2.7) 

2 r c r 

J = S + L. (2.8) 

Here the orbital angular momentum L is 

L = rxp + vxq = L N + L S o (2.9) 
with the Newtonian orbital angular momentum 

L N = fir xv (2.10) 

and the spin-orbit term 

Lso = f ((2 + »j)^[rx(rxS)]^[vx(vxS)]} . (2.11) 

Formally the spin terms arising in these expressions are of order eS. Since the spin is of order S ~ r 2 e 2 V, the terms 
containing the spin in the Lagrangian, the acceleration, the energy and the orbital angular momentum are of order 
e 3 / 2 V/c higher than the corresponding Newtonian terms. 
The spin precession equation |l| is: 

S = (4 + 377)^3 L N x S . (2.12) 

Hence the magnitude S of the spin is constant. Note that up to this point the description is valid up to (post) 2 - 
Newtonian order. 

Since we are interested in the leading spin terms of the gravitational radiation backreaction losses, we may replace 
Ln — > J by inserting terms of higher order in S. This is allowed, as we are neglecting terms of order eS 2 = e 5 in the 
precessional angular velocity, which are e 3 / 2 order higher then the spin terms that we have kept. (We shall return to 
the order-of-magnitude estimates at the end of this sec tion .) 



From the conservation of total angular momentum ( |2.8| ) and the (2.12) precession of S about J, the vector L is 
found to precess about J: 

L = -S = (4 + 3r7)-^JxL . (2.13) 



Combining ( 2.1 2| ) (with J instead of Ln) and (2.13) it follows that the angle k subtended by S and L is constant. As 



shown in Appendix A, the angle k is constant to order e 3 / 2 . Thus we have the picture of a paralelogram of constant 
sides L and S rigidly rotating about its diagonal J (Fig.l). 
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We choose the Cartesian coordinates r = {%, y, z} of the reduced-mass particle with origin at the center of mass 
such that J points along the z axis |l7f . In this approximation we can choose the conserved quantities characterizing 
the orbital motion as: the energy E (due to the time-independent nature of the Lagrangian), the magnitude L of L 
(the vector L undergoes a pure rotation in the order we are considering) and L$, the spin-oriented component of L 
(cf. Appendix A): 



E = L f 



L = 



(2.14) 



This description is close to the one in the Lense-Thirring approximation. The inclusion of finite mass does not destroy 
the basic features of the formalism developed in the test particle limit. In computing the expressions of the constants 
of motion, Ln • S was replaced by LgS, which is compatible with the order we are considering. Using polar coordinates, 
x = r sin 9 cos tp, y = r sin9 sirup, z = rcos9, the constants of the motion are: 



E = 



(IV 



Gmfi 



L 2 = l x 2 r\9 2 



^GLsS = ^. 2 + r 2 ( ^ +sill 2^ 2)] 
r c r I 

sin 9 tp ) — 4 - 1- — — ELso , 



Gm\x GLsS 
— + 7 i—r-r 



L, 



S 
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L COS K . 

From the first expression of the energy, v 2 can be expressed in terms of constants of motion and r: 



= 2 -E 



2Gr 



V- 



2GL S S 
c 2 fj,r 3 



(2.15) 
(2.16) 
(2.17) 

(2.18) 



From the second expression of the energy in ( 2.15 ) and from ( 2.16 ), the equation for the radial motion is found to 
decouple from the angular degrees of freedom: 



E 



Gm 



fi 2 r 2 



2-q 



EL S S 
c 2 mii 2 r 2 



2(2 + 77) 



GL S S 
c 2 fj,r 3 



(2.19) 



We conclude this section by considering some order-of-magnitude estimates in the approximation which has led to 
the picture of a rigidly rotating paralelogram. The neglected time-dependent terms in L 2 and L$ are (4 + 2>rj) J L ■ 
(L S o x S)/r 3 dt ps r 4 e 6 and (4 + 3ri)/2 JL- (S x L so) / r 3 Sdt « r 2 e 4 , respectively. These are of order e 3 / 2 and e smaller 



than the terms we keep. In the radial equation ( [2.19 ), as declared, we keep terms linear in the spin and drop terms 
of order e 3 / 2 smaller. 



III. PARAMETERIZATION OF THE ORBIT 



In this section, we employ Eq. (2.19) for determining the turning points and then finding convenient parametriza- 
tions of the orbit. We describe the procedure only in its outlines, as a detailed exposition was previously given for the 
Lense-Thirring case in I. The condition for turning poits r = yields a cubic equation. By looking for roots slightly 
different from the roots known for the nonspinning case, one unphysical root and the turning points r^ax are found: 
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Here Aq is the length of the Runge-Lenz vector to the zeroth order in the spin: 



(3.1) 



A Q = \ G 2 m 2 fj? + 



2EL 2 
A* 



(3.2) 



As in the Lense-Thirring case, we define the eccentric anomaly parametrization of the Kepler orbits, 
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c m/iAo 



cos£ , 



such that r™« is at £ = ir and £ = 0, respectively. The parameter derivative dr/dt, follows from ( |3.3[ ). Expressing 
1/f from ( 2.1 9| ) and linearizing in S 1 , the expression for dt/dt; is found: 



(3.3) 
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dt 

di 



1 dr 
7d£ 



fi 2 (Gmu - A Q cos£) ((2 + r))G 2 m 2 fj? - r)EL 2 )L s S cos £ 



iL 2 A {-2uE)- 



Integration of (3.4) from to 2n gives the orbital period: 

Gm/j, 3 



T = 2tt- 



(3.4) 



(3.5) 



Note that the expression for the period has the same functional form as in the Lense-Thirring case. 

The eccentric anomaly parametrization r = r(£) was useful in computing the orbital period, but it leads to unnec- 
essarily complicated expressions when instantaneous losses need to be averaged. For the latter purpose we introduce 
the true anomaly parametrization r — r(x) requiring that it has the following properties: 



(a) r(0) = 



d(cosx) 



and r(jr) = r, n 
-(70 + S^\)r 2 



(3.6) 
(3.7) 



where 70,71 are constants. Property (b) generalizes Kepler's second law for the area. The unique parametrization 
satisfying both (a) and (b) is: 



2(2 + n)GL s S A (2G 2 m 2 > 3 + EL 2 ) + Gm^(2G 2 m 2 fi 3 + 3EL 2 ) cos X 



ju(Gm/x + A cos x) 



c 2 L 2 A 

2r]ELsS Gmu 2 A + (G 2 m 2 u 3 



{Gran + A cos \) 2 
- EL 2 ) cosx 



c 2 mu 2 A{) 



{Gm^i + A cosx) 2 



In the same way as dt/d£ [Eq. (3.4 )] was obtained, dt/dx is found to have the form 
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dt 
dx 



1 dr 
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c 2 mL 4 



(2 + n)Gmn 2 (3Gmu + A cosx) - V EL ' 



(3.8) 



(3.9) 



The true anomaly parametrization (3.£) will be used for the integration over one period of all instantaneous losses. 
In the next section, all of these losses will be expressed in the form F = F( X )- They contain no other X dependent 
factor in the denominator than r 2+n , where n is positive integer. Time integration can be replaced by parameter 
integration; 



F(t)dt 



„, n dt , 
dx 



(3.10) 



where dt/dx is given by (3^). Here we encounter an advantage of property (b) of our parametrization: the x 
dependence of the denominator of the integrand is especially simple, r™. 



As in the Lense-Thirring case, the integrals are evaluated by computing the residues enclosed in the circle £ 



We find, as the second advantageous feature of the parametrization th at th ere is only one pole, at £ 
over one period is achieved by dividing the result by the period T [Eq.(3.5)]. 



0. Averaging 



IV. INSTANTANEOUS RADIATIVE LOSSES 



In this section, we obtain the instantaneous losses in the energy E, magnitude L and spin projection L$ of the 
orbital angular momentum. We then find suitable angular variables such that all losses can be rewritten in terms of 
the true anomaly parameter x- 

From the radiative power and the instantaneous total angular momentum losses of the binary system given by 
Kidder formulae (3. 24), (3. 28)], we keep the Newtonian and spin-orbit terms: 
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We want to compute the instantaneous radiative losses of the constants of motion. As we argued in the Introduction 
and will show in detail later in this section, the total angular momentum loss can be taken equal to the orbital angular 
momentum loss, dJ/dt = dh/dt. The loss in the magnitude of the orbital angular momentum is given by 

2L dL _ dL 2 _ djLjLj) _ 2L dLj 
dt dt dt dt 



The loss in is given as: 

dLs d 
~df ~ ~dt 



S\ _ dL S L dS f L S 

'Sj~~dt's + S'~dt~(s's 



dt) 



(4.4) 



Although the radiation-induced change in the spin S is e 1 / 2 order smaller than the spin-orbit part of the change in 
the orbital angular momentum, due to the vectors "L/S k e -1 ' 2 , the second and third terms are of comparable order 
with the first term. We will evaluate the terms depending on the structure of the spinning body later in this section. 
The first term is found simply by multiplying ( \L2[ ) by S/S. 

The next step is to find variables in which these radiative losses can be expressed solely in terms of the parameter 
X and constants of the motion. The polar coordinates (9, cj)) are not suitable for this purpose: 6 is not a monotonous 
function of time, thus the root of 6 2 can not be extracted unambiguously. 

We introduce new variables conveniently characterizing the separation vector r, the spin vector S and the orbital 
angular momentum vector L. The relative coordinates can be expressed in terms of the time dependent Euler angles 

in, $ as: 



x = r(cos $ cos ^ — cos ln sin $ sin^) , 
y = r(sin <I> cos ^ + cos in cos $ sin vp) , 



(4.5) 



z = r sin i n sin 4" 



Here in is the polar angle and $ — ir/2 the azimuthal angle of the Newtonian orbital momentum Ljv (Fig. 1). The 
variable ^! is the angle subtended by the momentary position vector and the node line (the intersection of the the 
plane orthogonal to the L^v with the plane orthogonal to the total angular momentum J). From the leading terms 
of the equations of motion (Cf. Appendix of I) it follows that \& is a monotonous function of t. Thus we extract the 
square root of the equation for ^E' 2 , choosing the positive root. 
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In the absence of spin, the orbital plane is at rest, and 'J is the usual polar angle, satisfying the area law. Thus the 
time derivatives of the angles are given by 



(4.6) 



We will not need the explicit expressions of the first-order (in S) terms \E , i, | i > i and l\ for parameterizing the instan- 
taneous losses. Rather it will be sufficient for us to find how they are interrelated. First, by computing the square 



L 2 N of the Newtonian orbital angular momentum in two different ways: 



L 2 N = (fir x v) 2 = L 2 - 2L 



N 



J SO 



where 2Ln • Lso are the spin terms from (2.16), the following relation is found: 

2Gmfi — rjEr 



*i = — cos tjv$i + cos K- 



c 2 m(ir 3 



Second, we compute the z component of Ln in two independent ways: 

(L N )z = (F x v) z = (L 2 - 2L N • L so )^ cosl n 



Linearizing the second equality of (O) in the spin, we obtain 

tan* 
$1 = H 

sm ijy 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



When the relations (4.8) and (4.10) are inserted in the instantaneous losses of energy, magnitude and spin projection 
of the orbital angular momentum, all terms with <3?i, l\ will cancel. 

The paralelogram spanned by the spin and orbital angular momentum is next described in terms of the new variables 
(Fig.l). We denote the polar and azimuthal angles of S and L by a, and i,0 + tt. During the first period when J 
lies on the z axis, K = I + a. Similarly as with L and Ln, the azimuthal and polar angles differ only by first order 
terms in the spin. As a result, no dependence on the azimuthal angles $ and (3 remains in the losses. 

As we mentioned in the Introduction, the spin is of e 1 / 2 order smaller than the orbital angular momentum. This 
implies that the side S of the paralelogram shrinks. As the angles i, ijv and a are present in the losses only in the 
spin-terms, we need only their zeroth order expressions: 



IN 



0, a = k . 



(4.11) 



The nonradiativc evolution of the angles i and tjv is discussed in Appendix B. 

From (2.18) and (2.19) v 2 and r 2 can be replaced by expressions containing solely r and constants. Making use of 
the foregoing observations, we obtain the instantaneous losses of energy, magnitude and spin component of the orbital 
angular momentum: 



dE 
~dt 



dL 
~dt 



O <_r III , 

8G 3 mLS cos n r 



2Gmfi 2 r + 11L 2 ) 



15c 7 fir s 
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8G 2 mL 
5c 5 /zr 5 
8G 2 Scosk 



(2fiEr 2 - 3L 2 



+ — ^^ {uGm^Er 3 + 3fir 2 (6EL 2 + G 2 mV) - llGm^tfr 
15c 7 fi 2 r< I 



-r) [2fi 2 E 2 r 4 + YlGm^Er 3 + 3[ir 2 (5EL 2 + G 2 m 2 ^) - 5Gm^ 2 L 2 r + 15L 4 ] | 



(4.12) 



dLs 
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5c 5 /ir 5 v 7 15c 7 /z 2 r 7 I 

47?[2//£ 2 r 4 + 12Gmfi a Er 3 + 3[ir 2 (hEL 2 + G 2 m 2 ^ 3 ) - 5Gm^ 2 L 2 r + 15L 4 ] 
+ sin 2 K^24Gm^ 3 Er 3 + 6^r 2 (6EL 2 - G 2 mV) + 72Gmfi 2 L 2 r - 90L 4 
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+rj[-Ay 2 E 2 r i - 24Gm^ 3 Er 3 + 6^r 2 (5EL 2 - G 2 m 2 fi 3 ) + 49Gm^ 2 L 2 r - 135L 4 ] \ 
fiLrr sin 2 k sin 2*| 36/xSr 2 + \2Gm[i 2 r - 18L 2 + rj (34/iSr 2 + \2Gm,[i 2 r - 15L 2 )} 
sin 2 KCOs2*|24Gm^ 3 i;r 3 + 6^r 2 (6.BL 2 + G 2 toV) + 6Gmfi 2 L 2 r + 18L 4 

+T][Afi 2 E 2 r 4 + 2AGm^L 3 Er 3 + 2[ir 2 {YiEL 2 + 3G 2 m 2 fi 3 ) + 3Gm^i 2 L 2 r + 15L 4 ] } j 



L dS\ /L S\ /S dS\ 



S di/ \S S J \S dt J 

Now we proceed to evaluate the two radiative spin-loss terms appearing in the Lg-loss. We need these terms 
to the Newtonian order and we compute them from the radiation-reaction potential. To lowest order, this is the 
Burke-Thorne potential Jjj||| 

V=~|f$W. (4-13) 

(5) 

where I- k is the fifth time derivative of the system's quadrupolc-moment tensor and yj are coordinates centered on 
the spinning body. The torque r is the integral of the backreaction force over the volume of the body: 

Ti = ^tiulf J P(y)y k y dv . (4.14) 

The torque induces the radiative change in the spin. The volume integral can be expressed with the tensor of inertia 
©jfc = / PiVjVk - S jk yiyi)dV: 

~dT = n= 5^ €iki y jk - ( } 

For an axisymmetrically spinning body the tensor of inertia is Ojk = diag(<3, 0, 0') in the system of principal axes 
with S on the z axis. The magnitude of the spin is S = 0'fi, where f2 is the angular velocity of the body. Thus 

1 dSi 2G , 6 , ( 5 ) 6 fi 

s^r = ^&- 1 ^ kI ^ SkSl (4 - 16) 

is perpendicular to the spin S. Here the quantity S = 0/0' — 1 is the deviation from sphericity and S the spin 
direction. When the body is rapidly rotating (V ^> e 1 / 2 ), it is centrifugally flattened (5 is nonnegligible) and the 
estimate of ACST 0j for dS/dt holds. For slow rotation, V rs e 1 / 2 , the deviation from sphericity is 6 is small. In both 
cases, th e app roximation above, equating dJ/dt with dL/dt is valid. 

From (4.16), the last term in the L5-I0SS vanishes, and the other term containing the radiative spin loss is 

I ■ f = - ^^ (1 - ^ ^ Sln2 * ( 12 ^ + 20Gm ^ r + 45L2 ) ^ 

+ 4L cos 2* (l8fj,Er 2 + 20Gmfi 2 r - 15L 2 ) | . 

The losses in the energy and in the magnitude of angular momentum contain constants of the motion and they de- 
pend only onr = f(x)- These expressions will take the suitable form for averaging, after inserting the parametrization 
(3.8). The loss dLs/dt depends also on the zeroth order expressions of the angle variable \& and of r, which are: 

f = t +x, f=^sinx- (4.18) 

With these expressions, the task of parametrizing is completed: the loss in L$ can also be expressed in terms of the 
true anomaly parameter \. 
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V. AVERAGED RADIATIVE LOSSES 



The ave raged losses are computed by the residue theorem, passing to the complex variable ( — e %x . The average of 
the term (4.17) vanishes, 



L dS\ 



(5.1) 



The averaged losses in the constants of motion are: 



(-) 
\dt I 



G 2 m( 2E ^) 3/2 {ugE 2 L 4 + 732 G 2 m 2 [i 3 EL 2 + 425G 4 mV) 
G 2 (-2E^ 2 Scoan 



(-) 

\dt I 



4G 2 m(-2£/i) 3 / 2 

5c 5 L 4 
G 2 (-2E^ 2 Scosk 
l5c 7 L 7 



520E 3 L 6 + 10740G 2 m V# 2 L 4 + 2A990G 4 m 4 fi 6 EL 2 + 12579G 6 mV 
+r 1 (256E 3 L 6 + 6660G 2 m 2 fi 3 E 2 L 4 + 16660G 4 m 4 ^ 6 ^i 2 + 8673G 6 toV; 
(UEL 2 + 15G 2 mV) 

1188£ 2 £ 4 + 6756G 2 mV££ 2 + 5345G 4 rnV 
■rj(772E 2 L 4 + U76G 2 m 2 ^ 3 EL 2 + 3665G 4 m 4 /^ 6 ; 



(5.2a) 



(5.2b) 



dLs 
dt 



4G 2 m(-2^) 3 / 2 cosK, „ - , G 2 (-2Ey) 3 ' 2 S 

-j£ 4 ^ EL + 15Gm ^ + 15 

| [ll88£ 2 L 4 + 6756G 2 mVi;£ 2 + 5345G 4 ?tiV + t](772E 2 L 4 + 4476G 2 mV.E£ 2 + 3665G 4 mV)] 



sin 2 k 



[3516E 2 L 4 + 17676G 2 m 2 /i 3 £:i 2 + 12975G 4 mV 



+r ? (2428E 2 L 4 + 122WG 2 m 2 (J?EL 2 + 9125G 4 mV)] 
sin 2 Kcos(2* ) [(26EL 2 + 33G 2 m 2 ^ 3 ) (12 EL 2 + 6G 2 mV) 

+r 1 (119EL 2 + 156G 2 m 2 (i 3 )(2EL 2 + G 2 mV)] X 



(5.2c) 



The expressions (5.2a) and ( 5.2b| ) agree with the available results of for the averaged power and < dL/dt >, 
despite the differences encountered in the higher-order terms of the definition of the e nerg y and angular momentum. 

The averaged loss in < dLg/dt > contains terms proportional to cos(25 , ). Ryan argues that in most cases 
such terms are averaged out by the Schwarzschild precession. 



VI. AVERAGED LOSSES IN TERMS OF ORBIT PARAMETERS 



As in I, we define the generalized semimajor axis and eccentricity imposing the conditions: 

r™ax = a(l ± e) . 



Solving the system and (6.1), one has: 

Gmfi 

a = 

2E 

8L S S 



1 



2EL S S 
c 2 mL 2 



(2 + ry) 



1 - e z = - 



2EL 2 
G 2 m 2 ii 3 



1 



2 „...._>..!, , 2r ' LsS ^EL 2 + 2G 2 m 2 f i 3 ) 



(EL 2 + G 2 m 2 \x 3 ) 



c 2 mL 4 c 2 mL 4 

From these relations one can express the constants of motion in terms of the orbit parameters 

Gmp, f | GScosk(2 + ??) 
2a ' 



L 2 = Gmn 2 a(l - e 2 )(l 



c 2 a 3 / 2 v /Gm(l -e 2 ) 

2GScosk 3 + e 2 + 2r? 
c 2 a 3 l 2 ^G~M\ -e 2 ) 1-e 2 



(6.1) 



(6.2) 



(6.3) 
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Inserting (3.3) in the equations (5.2a) we get the averaged radiation losses in terms of a, e and k: 
G 4 m 3 n 2 



(-) 
\dt I 



die 



15c 5 a 5 (l-e 2 ) 7 / 2 

G 9/2 TO 5/2 M 2 5c0SK 



(37e 4 + 292e 2 + 96) 

491e 6 + 5694e 4 + 6584e 2 + 1168 + 7?(355e 6 + 5316e 4 + 7248e 2 + 1200) 



dt 



30c 7 a 13 / 2 (l-e 2 ) 5 

(7e 2 + 8) 

549e 4 + 1428e 2 + 488 + r ? (403e 4 + 1366e 2 
G 4 m 2 n 2 S 



4^7/2^/2^2 ^ 

5c 5 a 7 /2(i - e 2 ) 2 
G 4 m 2 fi 2 S cos k 



15c 7 a 5 (l-e 2 ) 7 / 2 
4G 7 / 2 m 5 / 2 fi 2 cos « 



456) 



(6.4) 



(7e 2 + 8) + 



5c 5 a 7 / 2 (l-e 2 ) 2 v ; 30c 7 a 5 (l-e 2 ) 7 / 2 

| 2[549e 4 + 1428e 2 + 488 + 7 ? (403e 4 + 1366e 2 + 456)] 
- sin 2 K[l383e 4 + 4368e 2 + 1464 + 7 ? (1027e 4 + 3922e 2 + 1296)] 
+ sin 2 k cos(2* ) e 2 [l56e 2 + 240 + ?7(119e 2 + 193)] | . 



The losses of k, a and e follow by taking the time derivatives of (2.17) and of (3.2): 

G 7 /2 m 3/2 M 5 sinK 



/dn\ 
\~dt) 



ida\ 
\dt/ 



\dt/ 



30c 7 a 11 /2(i _ e 2)4 



|285e 4 + 1512e 2 + 488 + 77(221e 4 + 1190e 2 + 384) 
-e 2 cos(2* ) [I56e 2 + 240 + ?7(119e 2 + 193)] } 



G 7 / 2 m 3 / 2 fj,S cos k 



(6.5) 



15c 7 a 9 / 2 (l-e 2 ) 5 

G 3 TO 2 ^i 

15c 5 a 4 (l - e 2 ) 5 / 2 
G 7 / 2 m 3 / 2 /xScosK 
30c 7 a n / 2 (l-e 2 ) 4 



363e 6 + 3510e 4 + 7936e 2 + 2128 + ?7(291e 6 + 4224e 4 + 7924e 2 + 1680) 
e(121e 2 + 304) 

1313e 4 + 5592e 2 + 7032 + r;(1097e 4 + 6822e 2 + 6200) 



From Eq. (6.5), the angle k increases. Hence we find that the vectors L and S of a finite-mass binary system tend to 
antialign as in I. Substituting rj — for a test particle with negligible mass, we reproduce the Lense-Thirring losses 
of I. 

The evolution of the angles l and ln under radiation backreaction is determined by < dL/dt >. This computation, 
relegated to Appendix B, involves higher-order terms in S which are, however of (post) 3 / 2 -order. 



VII. CONCLUDING REMARKS 



The computation of finite-mass effects has proved to be much more difficult than the corresponding description of 
a test particle. We have been helped by several fortunate circumstances in overcoming the obstacles in computing 
the radiative losses up to the leading spin terms. First, in this approximation, one can still decouple a radial equation 
of motion. Furthermore, a generalization of the true anomaly parametrization x, found in I for the Lense-Thirring 
case, exists with the inclusion of the finite mass contributions. Although the dynamics of the system is much more 
complicated, we succeeded in describing all quantities entering in the radiative losses in terms of the parameter \- 

We were able to separately describe the relative evolution of the spin and orbital angular momentum vectors, up to 
first order in spin, under radiation backreaction. In the loss of the spin projection of the orbital angular momentum 
two terms, containing the radiative spin loss, appear. One of them is known to vanish (ACST). We computed the 
other term, and have shown that it averages out to nothing. We did not consider no-spin post-Newtonian effects 
here. In the equations of motion, we did not keep terms of higher order than e 3 / 2 , which, however, are comparable or 
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even larger than the radiative losses. They are needed for a complete characterization of the orbit, but they do not 
contribute to the radiative losses in the order to which we are computing them. 

In the test-particle (r) = 0) and constant-spin case our results agree with those of I. The Lense-Thirring description 
implies S ^> L. These pictures arise as different limits of a more general description in which both S and L are 
arbitrary. Such a description is subject to further investigation as it will go beyond considering only the linear terms 
in S. 
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APPENDIX A: NONRADIATIVE EVOLUTION OF L s AND k 



In this Appendix, we consider the nonradiative evolution of Ls and the angle k to linear order in the spin, and we 
show that these changes are smaller than order-e 3 / 2 . 



From the conservation of the total angular momentum J and the spin precession equation (2.12), the time derivative 
of Ls = L • S has the form: 



L, 



4 + 3ry G 2 fi 2 S 
2r 3 c 4 " ' 



T) ( S\2n / Sw S\/2 + J7 

+ 2kH) ]-{ r -s){ y 's){-T L + 



2Gm, 



(Al) 



The multiplicative factor of S needs to be evaluated o nly t o Newtonian order. We write the non-radiative change in 
Ls in terms of the radial and Euler-angle coordinates (4.5): 



Ls 



GS(4 + 37?) sin 2 k 



L(—Er)fxr 2 + 2Gm\i 2 r + L 2 r\) sin (2*) - L 2 rr\[ir cos(2*) 



After parameterizing by x, the averaged expression over one period is: 



G{-2E^f/ 2 SiiAl sin 2 n sin(2# ) 
16c 4 L 4 m 



(15ry 2 + 2677 ■ 



(A2) 



(A3) 



This expression is of order re 9 / 2 U, comparable with the leading radiative loss in L given by Peters and Mathews, 
corresponding to 2.5PN relative order. However we pursue our description to the 1.5PN relative order. If further 
terms were considered, e.g. e 5 / 2 , they would ge nera te additional terms in L, as well as in Ls- It is simple to check 
that these terms would be of the same order as (A3). Thus the change in Ls is negligible in our approximation. 
Similarly the change in cosk is negligible, as follows from (cosk)'= Lg/L. 



APPENDIX B: EVOLUTION OF 1 AND l n 



In this Appendix, we discuss the evolution, both in the absence of radiation and under radiation backreaction, of 
the angles 



L J 

cos I = — ■ — 



COSiAT = 



J 



(Bl) 



where J and Ln are the magnitude of the total angular mom entu m J and Ln, respectively. 
An expansion of the expressions on the right hand sides of (Bl) to 1.5PN order yields 



COS l = COS Iff = 1 



I.9/5 
— sin k 



cos K sin K 



S 



(B2) 



By this relation, both the radiative and the nonradiative changes of cos 1 and cos tjv are given in terms of the 
corresponding changes in k, S and L: 
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s 

S cos l — 8 cos l 7v = — [ — 



2 r 



. , . 5S 5L 
cos k o cos k + sin k — — 



, / 5S 5L\ 

[1 — 3 cos k ) o cos k + 3 cos k sin n | — — I 



(B3) 



The first and second terms on the right hand side are of 1PN order and of 1.5PN order, respectively. 

We have shown in Appendix A that the nonradiative evolution of the angle k is negligible in our 1.5PN approxi- 
mation. Because L and S are also constants, (cost,)' — cosln — under the nonradiative evolution. 



Finally we turn our attention to the radiative evolution of the angles l and ln . The spin term correction in (6.4) 
to the leading order radiative loss of L and the loss (6.5) of the angle k are of 1.5PN order, as they originate in the 



spin-orbit term Lso- Thus these terms do not contribute to the losses (B3). Similarly the loss in the magnitude of 
the spin S was seen to vanish by (4.16). Hence the only loss contributing to the radiative changes in i and ln is the 
leading term in the loss of L: 



dcost dcosiAr 



dt 



dt 



1 + 3 COS K- 



S\ 1 dL 



L L dt 



(B4) 



As the factors multiplying dL / dt are constant during a period of revolution, the averaged rates can be simply obtained 
by use of the leading terms in Eq. (5.2b) for < dL/dt >. 
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